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ABSTRACT 

Field  equations  of  nonlocal  elasticity  are  solved  to 

determine  the  state  of  stress  In  the  neighborhood  of 

a lln.  crack  in  an  elastic  plate  subject  to  a uniform 

shear  at  the  surface  of  the  crack  tip.  A fracture  criterion 

based  on  the  maximum  shear  stress  gives  the  critical 

value  of  the  applied  shea-;  for  which  the  crack  becomes  unstable. 

Cohesive  stress  necessary  to  break  the  atomic  bonds  is 

calculated  for  brittle  materials. 

1 . INTRODUCTION 

In  several  previous  papers,  [1]  - [3],  we  discussed  the  state  of 
stress  near  the  tip  of  a sharp  line  crack  in  an  elastic  plate  subject 
to  a uniform  tension  perpendicular  to  the  line  of  the  crack  at  infinity.  The 
solution  of  this  problem  was  obtained  within  the  framework  of  the  non- 
local elasticity  theory.  The  resulting  solution  did  not  contain  the 
stress  singularity  present  in  the  classical  elasticity  solution  and 
therefore  a natural  fracture  criterion  based  on  the  usual  maximum  stress 
hypothesis  could  be  established.  This  most  Interesting  outcome  could 
be  used  to  calculate  the  cohesive  stress  in  various  materials. 

^Thc  present  work  was  supported  by  the  Office  of  Naval  Research 
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T1k>  prc'si'iit  paper  deals  with  the  problem  of  a line  crack  In  an 
elastic  plate  where  the  crack  is  subject  to  a uniform  shear  load.  We 
employ  the  field  equations  of  nonlocal  elasticity  theory  to  formulate 
and  solve  this  problem.  Grat  Ify Inf'ly  the  resulting  solution  does  not 
contain  the  stress  singularity  at  the  crack  tip  and  therefore  a fracture 
criterion  based  ca  the  maximum  shear  .stress  hypothesis  can  be  used  to 
obtain  the  critical  value  of  the  applied  shear  for  which  the  line 
crack  begins  to  beconu.'  unstable.  If  the  concept  of  the  surface  energy 
is  used  It  is  possible  to  calculate  the  cohesive  stress  holding  the 
atomic  bonds  together.  For  steel  (with  no  dislocations)  we  give  an  estimate 
for  the  cohesive  stress.  In  section  2 we  give  a resum^  of  basic 
equations  of  the  linear  nonlocal  elasticity  theory.  In  section  3 the 
boundary  value  problem  is  formulated  and  the  general  solution  Is 
obtained.  In  section  4 the  solution  Is  given  for  the  duel  Integral 
equations  completing  the  solution  of  the  problem  of  line  crack  subject  to 
a shear  load.  Calculations  for  the  shear  stress  are  carried  out  on  a 
computer  and  results  are  di.scussed  in  section  5. 


2 . BASIC  EQUATIONS  OF  NONLOCAL  ELASTICITY 
The  basic  equations  of  linear,  homogeneous,  isotropic,  nonlocal  elastic 
sol'ds,  with  vanishing  body  and  Inertia  forces,  are  (cf.  |4,  5)) 


(2.2) 


kS, 


[A'  ( |x'-x|  )e^^(x')6|^jj^+2p’  ( |x’-x|  )ej^j^(x')  ]dv(x') 


'ki 


2 


(2.3) 
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wtiere  tlie  only  ctlfforence  from  classical  elasticity  is  in  the  stress 
constitutive  equations  (2.2)  which  states  that  the  stress  tj^j^(x),  at  a 
poi"t  X,  depends  on  strains  ej^^(x'),  at  aljL  points  of  the  body.  For 
homogeneous  and  isotropic  solids  the  nonlocal  elastic  moduli  X'(|x'-xl) 
and  p'(|x'-xl)  are  functions  of  the  distance  between  the  variable  point 
x'  and  the  fixed  point  x at  which  the  stress  is  to  be  evaluated.  The 
integral  in  (2.2)  is  over  the  volume  1/  of  the  body  enclosed  within 
the  surface  31^. 

Throughout  this  paper  we  employ  cartesian  coordinates  Xj^  with 
the  usual  convention  that  a free  index  takes  the  values  (1,  /,  3)  and 
repeated  indices  are  summed  over  the  range  (1,  2,  3).  Indices  following 
a comma  represent  partial  differentiation,  e.g. 


3uj^/3x^ 


In  our  previous  work  [4,  6,  7]  we  have  obtained  the  forms  of 
A'(|x'-x|)  and  u'(|x'-x|)  for  which  the  dispersion  curves  of  plane 
waves  coincide  within  the  entiie  Brlllouin  zone  with  those  obtained  in 
the  Born-Von  Kdrmdn  theory  of  lattice  dynamics.  Accordingly 


(2.4)  (A',p')  = (X,y)a(|x'-x|) 

a(lx'-x|)  ■ I (a-[x'-x| ) 
0 


x'-x| <a, 


where  a is  the  la*^^tlce  parameter,  X and  p are  classical  Lame  constants, 

and  a is  a normalization  constant  to  be  determined  from 

o 
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(2.5) 


J a(  I X '-X  I )dv(x ' ) “ 1. 
1/ 


While  this  sltnpie  and  elegant  result  Is  useful  in  many  calculations, 
it  is  not  the  only  one  that  approximates  the  dispersion  curves  of 
lattice  dynamics.  In  fact  a very  useful  one  is 


(2.6) 


•x|) 


a^exp[-(B/a)  (^c'j^-Xj^)  (x’j^-Xj^)  ] 


where  6 is  a constant.  For  the  two-dimensional  case  (2.6)  has  the 
.specific  form 


(2.7)  a(|x'-x|)  - ^ (6/i)^exp{-(B/a)^[  (x' j^-Xj^)^  + (x'2-x,)^]} 


Employing  (2.4),  in  (2.2)  we  write 


(2.8)  t^^  - I a(|x'-xi)  aj^j^(x')dv(x’) 


where 


(2.9) 


is  the  classical  Hooke's  law.  Substituting  (2.8)  into  (2.1)  and  using 
the  Green-Gauss  theorem  we  obtain 


f a(|x'-x|)0|^^  j^(x')dv(x')  - I a(|x'-x|  )Oj^,(x')daj^(x')  =0 

(/ dl  ' ' ''  ' 


(2.10) 


Hero  the  surface  Integral  may  be  dropped  If  the  only  surface  of  the 
body  Is  at  Infinity. 


3 . CRACK  UNDl'R  SHEAR 

Consider  a pl.ite  in  (x^Hx,  ^2=7)  - plane  weakened  by  a lino  crack 
of  length  2i  along  the  x-axls.  The  plate  is  subjected  to  a constant 
shear  stress  along  the  surfaces  of  the  crack,  Fig.  1.  For  the 
plane  strain  problem  (2.10)  takes  the  form 


(3.1) 


I a( jx'-x 


a(  Ix’-xl ) [o2£(>^' .0)  )dx'  = 0 


where  the  Integral  with  a slash  is  over  the  two-dimensional  infinite 

space  excluding  the  crack  line  (|x|<S.,  y“0) . A bold-face  bracket  indicates  a 

Jump  at  the  crack  line. 

When  an  Incision  is  made  in  an  unde.'ormed  body,  the  body  will 
in  general  be  deformed  and  stressed  because  of  the  long-range  inter- 
atomic attractions.  Thus  if  we  are  to  treat  the  problem  of  a plate 
with  crack,  undeformed  and  unstressed  in  the  natural  state,  we  must 
consider  that  after  an  incision  is  made  the  crack  is  not  opened,  l.e. 
the  boundary  conditions  are  to  be  applied  to  the  plate  in  the  natural 
state. 

I'nder  the  applied  uniform  shear  load  on  the  unopened  surfaces  of 
the  crack  the  displacement  field  possess  the  following  symmetry  regulations 


(3.2) 


u(x,-y) 


-u(x,y) 


v(x,-y)  »■  v(x,y) 
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Employlng  this  in  (2.9)  we  see  that 
(3.3)  [O2^(x,0)l  = 0 , |x|>£ 


Hence  the  limits  in  the  second  integral  of  (3.1)  may  be 

replaced  by  (-”,®) . 

The  Fourier  transform  of  (3.1)  with  respect  to  x'  gives 


(3. A) 


00 

* 

a(C. 

-•00 


y'-y|)[-iC0jj^(^,y')  + ^ ' 


-a(5,|y|)t  ] - 0 


where  a superposed  bar  indicates  the  Fourier  transform,  e.g., 


f (S,y) 


•00 

f (x,y)exp(lf,x)dC 

boo 


If  we  take  the  Fourier  transform  of  (3. A)  with  respect  to  y,  and  solve 
for  the  fact^  of  a in  the  integrand  of  (3. A),  upon  inversion  we  obtain 


-lC0j^(f;,y) 


[ O2j^(x,0) 


since  the  left  hand  side  is  defined  for  all  y except  y=0.  Thus  we 
obtain 


.1,;  + dSLt 

*•111  dy 


t 


(3.5) 


I “ 1,2 
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Substltutlng  this  Into  (3.4)  we  have 

(3.6)  [ O2j(x,0)  1=0  . It  = 1,  2 

Hence  we  have  shown  that  the  general  solution  of  (3.4)  is  the 
same  as  that  of  the  system  (3.5)  and  (3.6). 

The  jump  condition  (3.6)  for  lxl>£,  is  satisfied  identically.  Kor 
|x|<'£  we  also  have  (o^j^(x,0)  ] = 0,  on  account  of  (3.2).  Since 
022^^‘~y^  = -022^^’^^  also  see  that  t22(x*0)  “ 0 foi*  |x|<£.  Thus 
the  normal  stress  condition  on  the  crack  surface  is  satisfied  identically. 
Considering  also  the  continuity  requirement  of  the  displacement  field 
satisfying  (3.3)  we  find  that  the  boundary  conditions  at  x=0  are: 

a (x,0)  =0  , t (x,0)  = T , |xl<£ 

(3.7)  yy  yx  ’ O ’ ' > 

Oyy(x,0)  = 0 , u(x,0)  = 0 , |x|>£ 

In  addition  we  must  have 


(3.8)  u=v=0  , asy-^® 

Consequently  we  must  obtain  the  solution  of  (3.5)  subject  to  (3.7)  and 

(3.8) .^ 


^Even  though  some  authors  feel  that  v(x,0)“0,  |xl>£  should  also  be 
satisfied  for  this  (so-called  Mode  II)  prt)blem  (cf.  [8]),  the  results 
based  on  the  boundary  conditions  (3.7)  and  (3.8)  are  accepted  by 
workers  on  the  theory  of  fracture.  The  physical  reasoning  indicates 
that  constant  shear  load  if  not  balanced  by  an  opposing  couple,  should 
give  a rotation  to  the  whole  body  so  that  v(x,0)=0  for  |x|>£  appears 
to  be  in  contradiction  with  the  expected  displacement  field. 
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Kquatlons  (3.5)  are  nont>  other  than  the  Fourier  transforms  oi  the 
Navler's  equationgin  two  dimensions,  namely 


(3.9) 


- (X+2y)^^u-iC(X+M)v,y  = 0 , 


-if;(X+u)u,y+(A+  2u)v,^^-f;^liv  ” 0 


The  general  solution  of  this  set  (for  y>0  ),  satisfying  (3.8)  is. 


u = -(2tt) 


-1/2 


(3.10) 


V =■  (2ti) 


-1/2 


D 

[U|a(0  + d^ly-  ^^)B(t:)]exp(-|c|y- 


i[A(^)  + yB(C)]exp(-Uly-iSx)df. 


where  A(0  and  B(0  are  to  be  determined  from  the  boundary  condltons 
(3.7).  Using  (2.9)  we  calculate 


(3.11)  Oyy(C.y)  = 2u1[-U|A(C)  + ( ^ - U|y)  B(0)exp(-|c|y) 


According  to  (3.7)^  and  (3.7).^,  this  must  vanish  at  y=*0.  Hence 


(3.12)  B(0  * ^ U|  A(0 


Noting  that  A(-C)  ■ A(£i),  on  account  of  symmetry  u(x,y)  u(-x,y), 
the  displacement  field  may  be  put  Into  the  form 


(3.13) 


u(x,y)  - (h  h+iL 

TT  P 


v(x.y)  - (-)^^2  J&M 


A(0  “^iyj  e~^^cos(5x)dC , 


A(r t-  5y|  e ^^sln(5x)df, 

■ \T+r  ) 


15x)dC, 
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For  t'ne  ^hrough  (2.9)  and  ( 5.12)  w('  obtain  (y>0) 


" -2(2/it)^^^(X+p) 


A ( f, ) ( 2 r.  - f, y ) o ■ y s i n ( f;  X ) d r , 


0.14)  * -n,  (x,-y)  = -2^2.0)^/^  (>+,)) 


,oo 


yy 


yy 


A (?,  )f,^ye  in  (f;x)di; , 


= -2(2/7t)^/^(X+p) 
yx 


A (O  ^(l-f;y)e  ^’^cos (f;x)df; 


The  stress  field  according  to  (2.8),  is  then  given  by 


dy' 


Oj^jj(x' ,y ' ) (u(|x'-x|  , |y'-y|  ) - oi  ( | x '-x  | , | y '+ y | ) | dx 


(3.15) 


tyy(X.y) 


dy' 


'’yy^’^^y')  ta(|x'-x|  , |y'-y|  ) - o ( | x ' -x  [ , | y '+y  | ) ) dx 


t (x,y) 
yx  ■' 


dy' 


^yx^’^'*/')  [a(|x'“x|  , |y'-y| ) + a ( | x '-x  1 , | y '+y  | ) ]dx 


Substituting  for  a from  (2.7),  the  Integrations  may  be  performed  with 
respect  to  x'  and  y'  by  noting  the  integrals  (cf.,  [9|). 


1 


1 


2 

exp(-px’  ) 

4 

mm  UD 


sin  ^(x+x') 
cos  r,  (x+x') 


dx' 


(ii/p)  ^'^^cxp(-C^/4p) 


sin  f X 
cos  f,  X 


(3.16) 


ton 

exp(-py'^-Yy')dy'  - y(ii/p)  ' '^^exp(Y^/4p)  [ l-<J>(Y/2/i7)  1 , 

0 


y'exp(-py'^-Yy*)dy'  * ^ (^-)  ^^^^xp(Y^/4p)  [ 1-<Ky/2/p)  ] . 


Hence, 
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t . - -(2/ti)^ 


XX 


(2/Ti)"(X+ii) 


A(C)  sln(^x)  {e~^’y[2^+(5^/2p)(5-2py)] 


* I l-e^y[2^+(f  ^/2p)  (£+2py)  ] f 1 )dC  . 

^OO 

t =.  (2/n)^(A+p)  k^A(r,)aln(Cx)  { 1-<J ] 

yy  J„  2p  2/p 


(3.17) 


‘•yx  “ -(2/n)^(X+p) 


A(5)Ccr>s(Cx)  (e“^^[l+(C/2p)  (C-2py)  ] 


•U-<K^J^)1  + e^yil+(5/2p)(£;+2py)][l-^(^|^)] 


-2?(irp)  ^exp[-py^-(5^/4p)  ] }d£; 


where  p*(6/a),  and  't'(z)  is  the  error  function  defined  by 


't>(z)  » 2tt 


-i 


^ ""  2 
exp(-t  )dt 

o 


The  remaining  two  boundary  conditions  and  (3.7)^  may  now  be 

expressed  as 


(3.18) 


I i i 

r,k(4)  K(;,p)  cos(Cz)  dc=  -(Tt/2)"T 

Jo  ° 

r -i 

^ C(C)  cos(4z)  di;  “ 0 

0 


0<Z  <1 


Z>1 


where  we  have  introduced 
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x/n  1 Z , CJI  = c , MO  = r^c;(0, 


(3.19)  K(r,.p)  H [l+(f,2/2p£^)][l-<H^/2il^)]-(r,/Jl/W^exp(-';^/^pil^)  , 


= t^£^/2(A+p)  , p H (e/a)^ 


To  determine  the  unknown  function  A(^)  we  must  solve  the  dual  integral 
equations  (3.18)  for  C(c). 


4.  THE  SOLUTION  OF  DUAL  INTEGRAL  EQUATIONS 
By  introducing 

cos(tz)  = (iTi;z/2)  J_^(cz) 


where  J (z)  is  the  Bessel’s  function  of  order  we  write  the  system 

V 

(3.18)  in  the  form 


^C(?)  [l+k(eO  jJ  ,(^z)dc  = -T  z ^ , 0<z<l 

Jo  ° 

(4.1) 

»0O 

C(C)J_,  (cz)d?  -=0  , z > 1 

0 * 


The  kernel  function  k(  ^ 5)  is  given  by 


(4.2)  k(cO  = K(C,p)-l  = 2e^C^[l-<t>(E01-4>(e0-27i"^E^exp(-c^(;^) 


where 


(4.3) 


E = l/2Jlp^-  a/28£ 
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The  solution  of  the  dual  integral  equations  (4.1)  is  not  known.  However, 
It  is  possible  to  reduce  the  problem  to  the  solution  of  a Fredholm 
equation  (see  [10]  § 4.6). 


(4.4) 


h(x)  + 


,1 

o 


h(u)L(x,u)du  = -H^tx) 


for  the  function  h(x),  where 


(4.5)  L(x,u)  = (xu)‘ 


tk(et)J  (xt)J  (ut)dt 

0 0 


Once  (4.4)  is  solved  then  C(C)  is  calculated  by 


(4.6)  C(C)  - (20 


x^J^(Cx)h(x)dx 


We  observe  that  for  c=0  we  have  k=0,  and  the  dual  integral  equations 

(4.1)  reduce  to  those  obtained  in  the  classical  elasticity.  For  this 

case  from  (4.4)  we  have  h (x)  =»  -T  (Tix)^/2.and  (4.6)  gives  the  classical 

o 0 ’ 

result ; 


(4.7)  C (O  = -(n/2)^T  c”^Ji(0 
0 o i 


or 


(4.8)  A^(0  - -(tt/2)^T^Jj(^)0/5JI 


The  next  iteration  of  the  Integral  equation  (4.4)  with  the  use  of 
h^ (x)  gives 
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(^.9)  h^(x)  = -(Ti/2)^T^j;'^[r’ji(c)  - 11^,0)] 


where 


(4.10)  I(;,c)  H 


ri 


X JqCCx)!  (x,f.)dx 


(4.11) 


[u^(t)J^(0-tJ^(OJ^(t)]dt 


0 C -t 


-1. 


Jn  which  f is  the  Hankel  transform  of  t~  k(et) Jj^(t) , i.e., 


f(x,e)  = 


tf  (t ,e) J^(xt)dt , 


-1, 


f(t,e)  = t k(et)lj^(t) 


If  we  write 


x/e  - y,  ec  ■=  n 


in  (4.10)iWe  will  have 


I(C,0) 


2 

0 f(ey,e)dy 

0 


2* 

But  since  c f(ey,e)  is  the  Hankel  transform  of  f(n/0,0),  we  see  that 


lim  I(C,e)  * lim  f(4,e)  « lim  (C)k(EC) ] 

0-^0  0-*-O  c-H)  ^ 


-1/,- 


1 

i 


Hence  (4.10)  vanishes  as  k(e(;)  with  c-*-0,  l.e., 
(4.12)  n^(0  « -(Tr/2)H^r^J^(0  0(4) 


Thus  for  small  z the  (difference  between  ^^^(4)  and  0^(4)  will  be  of  order 
c.  Since  any  number  of  iterations  will  contribute  higher  order  terms 
in  c,  the  solution  for  C(4)  will  also  be  in  the  form  (4.12)  as  t:-*0. 

We  observe  that  in  general  6 is  in  the  neighborhood  of  1,  [6). 

The  ratio  of  the  atomic  distance  to  the  crack  length  a/29,  is,  however, 

extremely  small  for  even  microscopic  cracks.  Thus  e is  very  small 

and  k(e4)  is  in  general  negligible  as  compared  to  unity  for  all  values 

of  g4.  As  can  be  seen  from  Fig.  2,  k(e4)  go^s  uniformly  from  0 to 

-1  as  E4  varies  between  0 and  ®.  Thus  we  expect  that  the  solution  of  (4.1) 

for  C(4)  will  be  almost  the  same  as  the  classical  solution  C (4) 

0 

corresponding  to  c=0.  Some  differences  are  of  course  expected 
for  crack  lengths  close  to  atomic  distances. 

Substituting  A^(0  for  A((^)  in  (3.13)  and  (3.17),  we  obtain  the 
displacement  and  stress  fields.  It  la  of  Interest  to  calculate  the 
shear  stress  along  the  x-axls.  This  is  given  by 


(4.13) 


t /t 
yx  o 


.30 

[ l+k(c4)  ] Jj^(4)cos(4z)d4 

0 


Z r.  X/I 


This  integral  for  0<z<l  gives  ■ 1 since  in  this  Interval  the 

Integral  converges  even  for  e«0.  For  z>l  the  Integral  again  converges 
for  all  E>0  and  it  is  permissible  to  Ignore  k(e4)  as  compared  to  unity. 
However  for  z*l  this  is  no  longer  the  case  and  we  cannot  ignore  k(e4) 
as  compared  to  unity.  To  see  this  we  write  t^^^  at  z“l  as: 


I 


< 
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* tj 


where 


h - 


Jj^(n)  [l-ifCen) Icosndn, 


0 

.00 


2t 


.2  2 


e n Jj^(n)  ( l-<t(en)  ] cosndn , 


-2t  71 
0 


-i 


2 2 

enJj^(n)exp(-e  n )cosndn 

0 


Since  l-4>(en)iO  and  T iO  it  is  clear  that 

0 


I r i 

< (t^/  e)  (l“«l>(y))dy  - f /r.  e , 


4 ' > 


2 i 

y (l-'j'(y))dy  “ , 


t J s (2t  /tt^g) 
J o 


2 i 

y exp(-y  )dy  = t /it  c 

o 


Hence 


< (^3-  e’^ 


Thus  we  see  that  t (x,£)  has  finite  value  for  c^O.  In  the  evaluation 

y ^ 

of  the  stress  field  near  therefore  we  cannot  ignore  the  function 

k(Et)  in  (4.13). 
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5 . NUMERICAL  ANALYSIS  AND  DISCUSSION 

Calculations  of  the  shear  stress  t^^,  given  by  (4.13)  along  the 
crack  line,  were  carried  out  on  computer.  The  results 
are  plotted  for  e-1/20,  1/50,  1/100  and  1/200  in  Figures  3 to  6.  For 
a crack  length  of  20  atomic  distance  (e-1/20)  the  results  are  not  very 
good.  However  for  a cract7§l^l00  atomic  distances  (Fig.  5)  we  can  see 
that  the  shear  stress  boundary  condition  ty^(x,0)«T^  for  |x|<£  is 
satisfied  in  a strong  approximate  sense.  The  relative  error  is  less  than 
li  °/o.  Hence  we  conclude  that  the  use  of  the  classical  A^(C)  given 
by  (4.8)  gives  satisfactory  results  for  crack  lengths  greater  than  100 
atomic  distances. 

The  stress  concentration  occurs  at  the  crack  tip, and  it  is  given 
by 

(5.1)  ty^(£.,0)/r^  ■ c//T  , G 5 8a/2«, 
where  c converges  to  about  -0.30,  l.e, 

(5.2)  c - -o;30 

We  now  make  the  following  significant  observations: 

(1)  The  maximum  shear  stress  occurs  at  the  crack  tip,  and  it  is 
f inite  ( eq . 5.1). 

(il)  The  shear  stress  at  the  crack  tip  becomes  infinite  as  the 

atomic  distance  a-+0.  This  is  the  classical  continuum  limit 


of  square  root  singularity. 
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(iii)  When  t (£,0)  •=  t ( = cohesive  shear  stress) , f racture  will 
yx  c — 

occur.  In  this  case 


(5.3)  - G 

0 c 


where 


(5.4)  G - (Ba/2c^)tJ 
c c 


Equation  (5.3)  Is  non  other  than  the  expression  of  the  Griffith  criterion 
for  brittle  fracture.  Note  that  we  have  arrived  at  this  criterion  via 
the  maximum  shear  stress  hypothesis.  The  present  criterion  of  fracture 
not  only  unifies  the  fracture  mechanics  at  the  macroscopic  and  microscopic 
scales, but  also  employs  the  natural  concept  of  bond  failure  in  the 
atomic  scale. 

(Iv)  The  cohesive  stress  may  be  estimated  1;'  one  employs  the 
Griffith  definition  of  surface  energy  y and  writes 


(5.5) 


a = K y 


where 

(5.6)  K - 8c^u/ttB(1-v) 
c 

Employing  the  values  of  y and  the  elastic  constants  for  steel 

Y - 1975  CGS  , P - 6.92x10^^  CGS 

V - 0.291  , a - 2.48  A°  , B - 1.65 


-LH- 


we  find  tliat 


t “ 1.04x10^^  CGS  , t /ii  - 0..i  > 
c c 


This  result  is  in  the  right  range  and  well  accepted  by  metallurgists 

based  on  other  considerations.  For  example  Kelly  [11]  gives 

t - 6.0x10^^  CGS  and  t /g  - 0.11. 
c c 


Acknowledgement ! The  author  is  indebted  to  Dr.  L.  Hajdo  for  carrying  out 
the  computations. 


The  altemiation  constant  0=*1.65,  used  In  these  calculations,  makes 
the  Fourier  transform  of  « coincide  with  the  dispersion  curve  of 
elastic  waves  in  Botn-Von  Karman  model  of  lattice  dynamics  (cf,  (b)). 
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